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Abstract 

Network virtualization is an important concept to overcome the ossifica- 
tion of today's Internet as it facilitates innovation also in the network core 
and as it promises a more efficient use of the given resources and infras- 
tructure. Virtual networks (VNets) provide an abstraction of the physical 
network: multiple VNets may cohabit the same physical network, but can be 
based on completely different protocol stacks (also beyond IP). One of the 
main challenges in network virtualization is the efficient admission control 
and embedding of VNets. The demand for virtual networks (e.g., for a video 
conference) can be hard to predict, and once the request is accepted, the spec- 
ification / QoS guarantees must be ensured throughout the VNet's lifetime. 
This requires an admission control algorithm which only selects high-benefit 
VNets in times of scarce resources, and an embedding algorithm which re- 
alizes the VNet in such a way that the likelihood that future requests can be 
embedded as well is maximized. 

This article describes a generic algorithm for the online VNet embedding 
problem which does not rely on any knowledge of the future VNet requests 
but whose performance is competitive to an optimal offline algorithm that has 
complete knowledge of the request sequence in advance: the so-called com- 
petitive ratio is, loosely speaking, logarithmic in the sum of the resources. 
Our algorithm is generic in the sense that it supports multiple traffic mod- 
els, multiple routing models, and even allows for nonuniform benefits and 
durations of VNet requests. 



1 Introduction 



Virtualization is an attractive design principle as it abstracts heterogeneous re- 
sources and as it allows for resource sharing. Over the last years, end-system 



virtualization (e.g., Xen or VMware) revamped the server business, and we wit- 
ness a trend towards link-virtualization: router vendors such as Cisco and Juniper 
offer router virtualization, and Multiprotocol Label Switching (MPLS) solutions 
and Virtual Private Networks (VPNs) are widely deployed. Also split architectures 
like OpenFlow receive a lot of attention as they open new possibilities to virtualize 
links. 

Network virtualization 1 15 ] goes one step further and envisions a world where 
multiple virtual networks (VNets) — which can be based on different networking 
protocols — cohabit the same physical network (the so-called substrate network). 
VNet requests are issued to a network provider and can have different specifica- 
tions, in terms of Quality-of-Service (QoS) requirements, supported traffic and 
routing models, duration, and so on. The goal of the provider is then to decide 
whether to accept the request and at what price (admission control), and subse- 
quently to realize (or embed) the VNet such that its specification is met while 
minimal resources are used — in order to be able to accept future requests. 

Virtual networks have appealing properties, for instance, (1) they allow to inno- 
vate the Internet by making the network core "programmable" and by facilitating 
service-tailored networks which are optimized for the specific application (e.g., 
content distribution requires different technologies and QoS guarantees than, live 
streaming, gaming, or online social networking); (2) the given resources can be 
(re-)used more efficiently, which saves cost at the provider side; (3) start-up com- 
panies can experiment with new protocols and services without investing in an own 
and expensive infrastructure; among many more. 

Due to the flexibility offered by network virtualization, the demand for virtual 
networks can be hard to predict — both in terms of arrival times and VNet dura- 
tions. For example, a VNet may be requested at short notice for a video conference 
between different stakeholders of an international project. It is hence mandatory 
that this VNet be realized quickly (i.e., the admission and embedding algorithms 
must have low time complexities) and that sufficient resources are reserved for this 
conference (to ensure the QoS spec). 

This article attends to the question of how to handle VNets arriving one-by-one 
in an online fashion [9 ]: Each request either needs to be embedded or rejected, and 
the online setting means that the decision (embed or reject) must be taken without 
any information about future requests; the decision cannot be changed later (no 
preemption). 

The goal is to maximize the overall profit, i.e., the sum of the benefits of the 
embedded VNets. We use competitive analysis for measuring the quality of our 
online algorithm. The competitive ratio of an online algorithm is a if, for every se- 
quence of requests a, the benefit obtained by the algorithm is at least an a fraction 
of the optimal offline benefit, that is, the benefit obtainable by an algorithm with 
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complete knowledge of the request sequence a in advance. 
1.1 VNet Specification and Service Models 

There are many service models for VNets [26], and we seek to devise generic 
algorithms applicable to a wide range of models. The two main aspects of a service 
model concern the modeling of traffic and the modeling of routing. 

Traffic We briefly outline and compare three models for allowable traffic. (1) In 
the customer-pipe model, a request for a VNet includes a traffic matrix that spec- 
ifies the required bandwidth between every pair of terminals of the VNet. (2) In 
the hose model IPT61I2T1 . each terminal v is assigned a maximum ingress bandwidth 
bin{v) > 1 and a maximum egress bandwidth b out (v) > 1. Any traffic matrix that 
is consistent with the ingress/egress values must be served, (3) Finally, we propose 
an aggregate ingress model, in which the set of allowed traffic patterns is specified 
by a single parameter X > 1. Any traffic in which the sum of ingress bandwidths 
is at most X must be served. 

The customer-pipe model sets detailed constraints on the VNet and enables 
efficient utilization of network resources as the substrate network has to support 
only a single traffic matrix per VNet. On the other hand, the hose model offers 
flexibility since the allowed traffic matrices constitute a polytope. Therefore, the 
VNet embedding must to take into account the "worst" allowable traffic patterns. 

Multicast sessions are not efficiently supported in the customer-pipe model 
and the hose model. In these models, a multicast session is translated into a set 
of unicasts from the ingress node to each of the egress nodes. Thus, the ingress 
bandwidth of a multicast is multiplied by the number of egress nodes ifTTl [T8l l22l 

In the aggregate ingress model, the set of allowable traffic patterns is wider, 
offers simpler specification, and more flexibility compared to the hose model. In 
addition, multicasting and broadcasting do not incur any penalty at all since inter- 
mediate nodes in the substrate network duplicate packets exiting via different links 
instead of having multiple duplicates input by the ingress node. For example, the 
following traffic patterns are allowed in the aggregate ingress model with parameter 
X: (i) a single multicast from one node with bandwidth X, and (ii) a set of multicast 
sessions with bandwidths f\, where Yli fi < 1- Hence, in the aggregate ingress 
model traffic may vary from a "heavy" multicast (e.g., software update to multi- 
ple branches) to a multi-party video-conference session in which every participant 
multicasts her video and receives all the videos from the other participants. 
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Routing We briefly outline three models for the allowed routing. (1) In tree rout- 
ing, the VNet is embedded as a Steiner tree in the substrate network that spans the 
terminals of the VNet. (2) In single path routing, the VNet is embedded as a union 
of paths between every pair of terminals. Each pair of terminals communicates 
along a single path. (3) In multipath routing, the VNet is embedded as a union 
of linear combinations of paths between terminals. Each pair of terminals u and 
v communicates along multiple paths. The traffic from node u to node v is split 
among these paths. The linear combination specifies how to split the traffic. 

In tree routing and single path routing, all the traffic between two terminals of 
the same VNet traverses the same single path. This simplifies routing and keeps the 
packets in order. In multipath routing, traffic between two terminals may be split 
between multiple paths. This complicates routing since a router needs to decide 
through which port a packet should be sent. In addition, routing tables are longer, 
and packets may arrive out of order. Finally, multicasting with multipath routing 
requires network coding HI. 

Packet Rate We consider link bandwidth as the main resource of a link. How- 
ever, throughput can also depend on processing power of the network nodes. Since 
a router needs to inspect each packet to determine its actions, the load incurred on 
a router is mainly influenced by the so-called packet rate, which we model as an 
additional parameter of a VNet request. If packets have uniform length, then the 
packet rate is a linear function of the bandwidth. 

Duration and Benefit The algorithms presented in this article can be competitive 
with respect to the total number of embedded VNets. However, our approach also 
supports a more general model where VNets have different benefits. Moreover, we 
can deal with VNets of finite durations. Therefore, in addition to the specification 
of the allowable traffic patterns, each request for a VNet has the following param- 
eters: (i) duration, i.e., the start and finish times of the request, and (ii) benefit, i.e., 
the revenue obtained if the request is served. 

1.2 Previous Work 

For an introduction and overview of network virtualization, the reader is referred 
to lfl5l . A description of a prototype network virtualization architecture (under 
development at Telekom Innovation Laboratories) appears in ll3~TTl . 

The virtual network embedding problem has already been studied in various 
settings, and it is well-known that many variants of the problem are computation- 
ally hard (see, e.g., J21[14]]). Optimal embeddings in the multi-path routing model 
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exist for all traffic models. In fact, in the customer pipe model, an optimal mul- 
tipath fractional embedding can be obtained by solving a multicommodity flow 
problem. In the hose model, an optimal reservation for multipath routing in the 
hose model is presented in [18]. This algorithm can also be modified to handle the 
aggregate ingress model. 

Offline algorithms for tree routing and single path routing lack edge capacities 
and have only edge flow costs. Namely, these algorithm approximate a min-cost 
embedding in the substrate graph without capacity constraints. In the hose model, 
constant approximation algorithms have been developed for tree routing ifTTl l22l 
l24l . In the special case that the sum of the ingresses equals the sum of the egresses, 
an optimal tree can be found efficiently, and the cost of an optimal tree is within 
a factor three of the cost of an optimal reservation for multipath routing ll25l[29l . 
Kumar et. al ll29l proved that in presence of edge capacity constraints, computing a 
tree routing in the hose model is NP-hard even in the case where bi n (v) = b ou t(v), 
for every node v. Moreover, they also showed approximating the optimal tree 
routing within a constant factor is NP-hard. 

Published online algorithms for VNet embeddings are scarce. In Il23ll30l . an 
online algorithm for the hose model with tree routing is presented. The algorithm 
uses a pruned BFS tree as an oracle. Edge costs are the ratio between the demand 
and the residual capacity. We remark that, even in the special case of online vir- 
tual circuits ("call admission"), using such linear edge costs lead to trivial linear 
competitive ratios [5 1. The rejection ratio of the algorithm is analyzed in 112311301 . 
but not the competitive ratio. The problem of embedding multicast requests in an 
online setting was studied in l28l . They used a heuristic oracle that computes a 
directed Steiner tree. The competitive ratio of the algorithm in [28 ] is not studied. 
In fact, much research has focused on heuristic approaches, e.g., [20] proposes 
heuristic methods for constructing different flavors of reconfiguration policies; and 
P4l proposes subdividing heuristics and adaptive optimization strategies to reduce 
node and link stress. In [4], an online algorithm is presented for the case of multiple 
multicast requests in which the terminals the requests arrive in an arbitrarily inter- 
leaved order. The competitive ratio of the online algorithm in [01 is O (log n ■ log d) , 
where n denotes the number of nodes in the substrate network and d denotes the 
diameter of the substrate network. 

Bansal et al. [8 ] presented a result on network mapping in cloud environments 
where the goal is to minimize congestion induced by the embedded workloads, i.e., 
to minimize the edge capacity augmentation w.r.t. a feasible optimal embedding. 
They consider two classes of workloads, namely depth-d trees and complete-graph 
workloads, and describe an online algorithm whose competitive ratio is logarithmic 
in the number of substrate resources, i.e., nodes and edges. Moreover, every node 
in the workload is mapped to a node in the substrate network, and every edge is 
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mapped to a single path in the substrate network between its (mapped) nodes. In 
contrast, we allow arbitrary workloads, a wide range of traffic models and routing 
models, specify the mapping of nodes, and focus is on revenue maximization. 

Circuit switching is a special case of VNet embeddings in which each VNet 
consists of two terminals and the routing is along a single path. Online algorithms 
for maximizing the revenue of circuit switching were presented in |5]. A general 
primal-dual setting for online packing and covering appears in ifTOl [T2l . In the 
context of circuit switching, the load of an edge e in a network is the ratio between 
the capacity reserved for the paths that traverses e, and the capacity of e. In the 
case of load (or congestion) minimization the online algorithm competes with the 
minimum augmentation [lOl HI. In the case of permanent routing, Aspnes 
et. al [3] designed an algorithm that augments the edge capacities by a factor of 
at most O(logn) w.r.t. a feasible optimal routing. Aspnes et. al also showed 
how to use approximated oracles to embed min-cost Steiner trees in the context of 
multicast virtual circuit routing. 

1.3 Our Contribution 

This article describes an algorithmic framework called GVOP (for general VNet 
online packing algorithm) for online embeddings of VNet requests. This frame- 
work allows us to decide in an online fashion whether the VNet should be admitted 
or not. For the embedding itself, an oracle is assumed which computes the embed- 
dings of VNets. While our framework yields fast algorithms, the embedding itself 
may be computationally hard and hence approximate oracles maybe be preferable 
in practice. We provide an overview of the state-of-the-art approximation algo- 
rithms for the realization of these oracles, and we prove that the competitive ratio 
is not increased much when approximate oracles are used in GVOP. Our framework 
follows the primal-dual online packing scheme by Buchbinder and Naor iflOl [121 
that provides an explanation of the algorithm of Awerbuch et al. 0. 

In our eyes, the main contribution of this article lies in the generality of the 
algorithm in terms of supported traffic and routing models. The GVOP algorithm 
is input VNet requests from multiple traffic models (i.e., customer-pipe, hose , or 
aggregate ingress models) and multiple routing models (i.e., multipath, single path, 
or tree routing). This implies that the network resources can be shared between 
requests of all types. 

We prove that the competitive ratio of our deterministic online algorithm is, in 
essence, logarithmic in the resources of the network. The algorithm comes in two 
flavors: (i) A bi-criteria algorithm that achieves a constant fraction of the optimal 
benefit while augmenting resources by a logarithmic factor. Each request in this 
version is either fully served or rejected, (ii) An online algorithm that achieves a 
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logarithmic competitive ratio without resource augmentation. However, this ver- 
sion may serve a fraction of a request, in which case the associated benefit is also 
the same fraction of the benefit of the request. However, if the allowed traffic pat- 
terns of a request consume at most a logarithmic fraction of every resource, then 
this version either rejects the request or fully embeds it. 

1.4 Article Organization 

The remainder of this article is organized as follows. We introduce the formal 
model and problem definition in Section 12 The main result is presented in Sec- 
tion [3] The algorithmic framework is described in Section HI Section [5] shows how 
to apply the framework to the VNet embedding problem and discusses the embed- 
ding oracles to be used in our framework under the different models. The article 
concludes with a short discussion in Section [6] 

2 Problem Definition 

We assume an undirected communication network G = (V, E) (called the physical 
network or the substrate network) where V represents the set of substrate nodes 
(or routers) and E represents the set of links. Namely, {u, v} G E for u, v 6 V 
denotes that u is connected to v by a communication link. Each edge e has a 
capacity c(e) > 1. In Section 15721 we will extend the model also to node capacities 
(processing power of a node, e.g., to take into account router loads). 

The online input is as follows. The operator (or provider) of the substrate 
network G receives a sequence of VNet requests a = {ri, • • •}■ Upon arrival of 
request rj, the operator must either reject rj or embed it. A request rj and the set 
of valid embeddings of rj depend on the service model. A VNet request rj has the 
following parameters: (1) A set Uj C V of terminals, i.e., the nodes of the VNet. 
(2) A set Trj of allowed traffic patterns between the terminals. For example, in the 
customer-pipe model, Trj consists of a single traffic matrix. In the hose model, 
Trj is a polytope of traffic matrices. (3) The routing model (multipath, single path, 
or tree). (4) The benefit bj of rj. This is the revenue if the request is fully served. 
(5) The duration Tj = [tj ,tj] of the request. Request rj arrives and starts at 

time tj and ends at time tj. 

The set of valid embeddings of a VNet request rj depends on the set Trj of 
allowed traffic patterns, the routing model, and the edge capacities, for example: 
(1) In the customer-pipe model with multipath routing, an embedding is a multi- 
commodity flow. (2) In the hose model with tree routing, a valid embedding is a set 
of edges with bandwidth reservations that induces a tree that spans the terminals. 
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The reserved bandwidth on each edge may not exceed the capacity of the edge. In 
addition, the traffic must be routable in the tree with the reserved bandwidth. 

If the allowed traffic patterns of a request rj consume at most a logarithmic 
fraction of every resource, then our algorithm either rejects the request or fully em- 
beds it. If a request consumes at least a logarithmic fraction of the resources, then 
the operator can accept and embed a fraction of a request. If an operator accepts an 
e-fraction of rj, then this means that it serves an e-fraction of every allowed traffic 
pattern. For example, in the customer-pipe model with a traffic matrix Tr, only the 
traffic matrix e • Tr is routed. The benefit received for embedding an e-fraction of 
rj is e ■ bj. The goal is to maximize the sum of the received benefits. 

3 The Main Result 

Consider an embedding of VNet requests. We can assign two values to the embed- 
ding: (1) The benefit, namely, the sum of the benefits of the embedded VNets. (2) 
The maximum congestion of a resource. The congestion of a resource is the ratio 
between the load of the resource and the capacity of a resource. For example, the 
load of an edge is the flow along the edge, and the usage of a node is the rate of 
the packets it must inspect. A bi-criteria competitive online packing algorithm is 
defined as follows. 

Definition 1. Let OPT denote an optimal offline fractional packing solution. An 
online packing algorithm Alg is (a, (3) -competitive if: ( i) For every input sequence 
a, the benefit of Alg (a) is at least 1/a times the benefit of OPT. (ii) For every 
input sequence a and for every resource e, the congestion incurred by Alg (a) is at 
most [3. 

The main result of this article is formulated in the following theorem. Consider a 
sequence of VNet requests {rj }j that consists of requests from one of the follow- 
ing types: (i) customer pipe model with multipath routing, (ii) hose model with 
multipath routing, or single path routing, or tree routing, or (iii) aggregate ingress 
model with multipath routing, or single path routing, or tree routing. 

Theorem 1. Let [3 = 0(log(\E\ ■ (max e c e ) • (max,,- bj))). For every sequence 
{rj}j of VNet requests, our GVOP algorithm is a (2, (3) -competitive online all-or- 
nothing VNet embedding algorithm. 

Note that the competitive ratio does not depend on the number of requests. 
The proof of Theorem Q] appears in Sections H] and [5] 
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4 A Framework for Online Embeddings 



Our embedding framework is an adaptation of the online primal-dual framework by 
Buchbinder and Naor lfTTl[T2l . We allow VNet requests to have finite durations and 
introduce approximate oracles which facilitate faster but approximate embeddings. 
In the following, our framework is described in detail. 

4.1 LP Formulation 

In order to devise primal-dual online algorithms, the VNet embedding problem 
needs to be formulated as a linear program (LP). Essentially, a linear program 
consists of two parts: a linear objective function (e.g., minimize the amount of 
resources used for the embedding), and a set of constraints (e.g., VNet placement 
constraints). As known from classic approximation theory, each linear program 
has a corresponding dual formulation. The primal LP is often referred to as the 
covering problem, whereas the dual is called the packing problem. In our online 
environment, we have to deal with a dynamic sequence of such linear programs, 
and our goal is to find good approximate solutions over time lfITl[T2l . 

In order to be consistent with related literature, we use the motivation and for- 
malism from the online circuit switching problem @ (with permanent requests). 
Let G = (V, E) denote a graph with edge capacities c e . Each request rj for a vir- 
tual circuit is characterized by the following parameters: (i) a source node aj € V 
and a destination destj G V, (ii) a bandwidth demand dj, (iii) a benefit bj. Upon 
arrival of a request rj, the algorithm either rejects it or fully serves it by reserving 
a bandwidth of dj along a path from aj to destj. We refer to such a solution as 
all-or-nothing. The algorithm may not change previous decisions. In particular, a 
rejected request may not be served later, and a served request may not be rerouted 
or stopped (even if a lucrative new request arrives). A solution must not violate 
edge capacities, namely, the sum of the bandwidths reserved along each edge e is 
at most c e . The algorithm competes with an optimal fractional solution that may 
partially serve a request using multiple paths. The optimal solution is offline, i.e., 
it is computed based on full information of all the requests. 

First, let us devise the linear programming formulation of the dual, i.e., of on- 
line packing. Again, to simplify reading, we use the terminology of the online 
circuit switching problem with durations. Let Aj denote the set of valid embed- 
dings of rj (e.g., Aj is the set of paths from aj to destj with flow dj). Define a dual 
variable yj£ € [0, 1] for every "satisfying flow" fjj € Aj. The variable yj i spec- 
ifies what fraction of the flow fj t £ is reserved for request rj. Note that obviously, 
an application of our framework does not require an explicit representation of the 
large sets Aj (see Section [5]). 
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minZj -1 + X T -C s.t. 
Zj.D J +X T -A j >Bj 
X, Za > 



(I) 



max Bj ■ Yj s.t. 
A 3 -Y 3 <C 

DyYj<l 

Yj>0 



(II) 



Figure 1 : (I) The primal covering LP. (II) The dual packing LP. 



Online packing is a sequence of linear programs. Upon arrival of request rj, 
the variables y 3 ^ corresponding to the "flows" fjj € Aj are introduced. Let Yj de- 
note the column vector of dual variables introduced so far (for requests n, . . . , rj). 
Let Bj denote the benefits column vector (&i,i, . . . , &i,|Ai|> •••■> bj,i, ■ ■ ■ , &j,|A|) r > 
where W, k : b^i = 6j & for every i, hence we abbreviate and refer to ba simply as 
bj. Let C denote the "capacity" column vector (c\, . . . , cn) T , where N denotes the 
number of "edges" (or resources in the general case). The matrix Aj defines the 
"capacity" constraints and has dimensionality N x Yli<j An entry (Aj) e ^ i ^ 
equals the flow along the "edge" e in the "flow" f^g. For example, in the case of 
circuit switching, the flow along an edge e by f^i is di if e is in the flow path, 
and zero otherwise. In the general case, we require that every "flow" fj £ incurs a 
positive "flow" on at least one "edge" e. Thus, every column of Aj is nonzero. The 
matrix Aj + \ is an augmentation of the matrix Aj, i.e., |A J+ i[ columns are added 
to Aj to obtain A/+i. Let Dj denote a 0-1 matrix of dimensionality j x Yl<i<j I 
The matrix Dj is a block matrix in which {D^ ui^ = 1 if i = i', and zero 
otherwise. Thus, Dj + \ is an augmentation of Dj\ in the first j rows, zeros are 
added in the new |A J+ i| columns, and, in row j + 1, there are zeros in the first 
J2i<j I Aj| columns, and ones in the last |A J+ i| columns. The matrix Dj defines 
the "demand" constraints. The packing linear program (called the dual LP) and the 
corresponding primal covering LP are listed in Figure Q] The covering LP has two 
variable vectors X and Zj. The vector X has a component x e for each "edge" e. 
This vector should be interpreted as the cost vector of the resources. The variable 
Zj has a component Zi for every request rj where i < j. 



4.2 Generic Algorithm 

This section presents our online algorithm GVOP to solve the dynamic linear pro- 
grams of Figured] The formal listing appears in Algorithm Q] 

We assume that all the variables X, Z, Y, (i.e, primal and dual) are initialized to 
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zero (using lazy initialization). Since the matrix Aj + \ is an augmentation of Aj, we 
abbreviate and refer to Aj simply as A. Let coIqa (A) denote the column of A (in 
fact, Aj) that corresponds to the dual variable yj^. Let j(J, £) = X T ■ co\rjMA), 
where the values of X T are with respect to the end of the processing of request 
r 3 -_i. It is useful to interpret £) as the X-cost of the "flow" fji for request j. 
Let w(j,£) = • colji(A), namely, w(j,£) is the sum of the entries in column 
(j, £) of A. Since every column of A is nonzero, it follows that w(j, £) > (and 
we may divide by it). 

Algorithm 1 The General all-or-nothing VNet Packing Online Algorithm (GVOP). 
Upon arrival of request rj : 

1. fj t i ^— argmin{7(j, £) : fjj G Aj} (oracle procedure) 

2. If 7(7, £) < bj, then accept rf. 

(a) yj,e <- 1. 

(b) z^bj-^jj). 

(c) For each row e do: 

x e ^x e ■ 2 A ^U,')/ C « + } nx ■ (2 A *,U,t)l*> - 1). 
w(j,£) 

3. Else reject rj (note that X, Zj,Yj have not been changed.) 



Definition 2. Let Y* denote an optimal offline fractional solution. A solution Y > 
is (a, (3) -competitive if: (i) For every j, Bj -Yj>^- Bj ■ Y*. (ii) For every j, 
Aj -Yj <[3-C and Dj • Yj < 1. 

The following theorem can be proved employing the techniques of ifTTI . 

Theorem 2. Assume that: (i) for every row e of A, max.j^ A e u^ < c e , (ii) for 
every row e of A, min^ ^4 e ,(j/) £ {0} U [1> °°)> (Hi) for every column (i, £) of A, 
w(i,£) > 0, and(iv) min^ bj > 1. Let [3 = log 2 (l+3-(max,- ^ w(j, £))-(maxj &.,■)). 
r/jg GVOP algorithm is a (2, (5)-competitive online all-or-nothing VNet packing 
algorithm. 

Proof. Let us denote by Primal j (respectively, Dualj) the change in the primal 
(respectively, dual) cost function when processing request j. 

We show that Primal j < 2 • Dualj for every j. We prove that GVOP produces 
feasible primal solutions throughout its execution. Initially, the primal and the 
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(?") 

dual solutions are 0, and the claim holds. Let x e denote the value of the primal 
variable x e when rj is processed, i.e., after Step [2c] and before the execution of 

Step [2c] upon the arrival of request Tj + \, in particular xf^ = for every e. If 
rj is rejected then Primal j = Dualj = and the claim holds. Then for each 

accepted request rj, Dualj = bj and Primal j = Yl e eE(j e)( x ^ ~ x e ^) ■ c e + Zj, 
where E{j,t) = {e 6 {1,...,N} : A e rjg\ / 0}. Step l2cl increases the cost 
X T ■ C = x e ■ c e as follows: 



£ {xP- X V- l ~>).c = x *~ 1} ' (2 Ae '^ /ee - 1) + -J^ ■ (2 A °^ C * - 1) 

e£E(j,t) e€E(j,e) L 



< 



w(j,£) 

D (^ 1) + ^)'^ )/Ce - 1 )- 

£ (»?- 1) +^)^ttD=7^ + l. 



e£E(J,l) 

where the third inequality holds since maxj^ A e tj^ < c e . Hence after Step [2b] 
Prima/, < 7 (j, + 1 + (6j - 7(i, £)) = 1 + 6,- < 2 ■ 6,-, 

where the last inequality holds since minj bj > 1. Since Dwa/j = 6j it follows 
that Primal j < 2 • Dualj. After dealing with each request, the primal variables 
{x e } e U {zi}i constitute a feasible primal solution. Using weak duality and since 
Primal j < 2 ■ Dualj, it follows that: Bj-Y* < X T ■ C + Zj ■ 1 < 2 • flj • YJ- 
which proves 2-competitiveness. 

We now prove /3-feasibility of the dual solution, i.e., for every j, Aj-Yj < f3 ■ C 
and Dj ■ Yj < 1. First we prove the following lemma. Let row e (A) denote the eth 
row of A. 

Lemma 1. For every j > 0, 

X U) > 1 . ( 2 ™vv e (A J ).K,/ Ce _ y _ 

(max i}i w{i,£)) 

Proof. The proof is by induction on j. 

Base j = 0: Since the variables are initialized to zero the lemma follows. 
Step : The update rule in Step [2c] is 

X e <r- X e ■ 2 A ^^) /Ce + - -r • (2^^)/^ - 1) . 

Plugging the induction hypothesis in the update rule implies: 
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xi j) = x^-^ -2 A "^^ /c ' +— L— . (2 A ^> /ee - 1) 

w(j,£) 

1 



(max^ «;(»,*)) ^ 1J 1 J 



(maxj^ w(i,£)) (max it ew(i,£)) 

I 2 row c (A 3 ).y,/ Ce _ 



(maxi,£ tu(i, £)) (maxj^ w(i,l)) 

The lemma follows. □ 

Since for every row e of A, vaxsij^A e u^\ G {0} U [1, oo), it follows that in 
Step[2c]it holds that for every e such that ^4 e ,0/) 7^ 0, 

< 6,- • 2 Ae >(^> /Ce H f2 Ae '0- f ) /Ce - 1) 

e J w(j,e) { 

Since for every row e of A, (maxj^ A e up,) < c e , min^ A e uQ G {0} U [1, oo), 
and since for every column (i, £) of ^4, it?(«, £) > 0, and (mirij &,) > 1, it follows 
that, 

4 i} < 2 • &j + 1 < 3 • bj . 

Lemma [T] implies that: 

• (2™^ A i)- Y i/ c e - 1) < Xe < 3 ■ bj < 3 • (max 6;) . 



(maxj i£ w(i,£)) 
Hence, 



row e (A,) • Yj < log 2 [l + 3 • (max io (i, £)) • (max&j)] • c e , 

for every j, as required. □ 

Remark 1. 77je assumption in Theorem\2\that maXj^Ag < c e for every row 
e means that the requests are feasible, i.e., do not overload any resource. In our 
modeling, if rj is infeasible, then r,j is rejected upfront ( technically, Aj = 0). In- 
feasible requests can be scaled to reduce the loads so that the scaled request is 
feasible. This means that a scaled request is only partially served. In fact, multiple 
copies of the scaled request may be input (see ^2/ for a fractional splitting of re- 
quests). In addition, in some applications, the oracle procedure is an approximate 
bi-criteria algorithm, i.e., it finds an embedding that violates capacity constraints. 
In such a case, we can scale the request to obtain feasibility. 
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If a solution Y is {a, (3) -competitive, then Y//3 is a ■ /3-competitive. Thus, we 
conclude with the following corollary. 

Corollary 1. The GVOP algorithm computes a solution Y such that Y//3 is a 
fractional 0(f3) -competitive solution. 

Consider the case that the capacities are larger than the demands by a logarith- 
mic factor, namely, min e c e //3 > max. j^g A e un for every row e of A. In this case, 
we can obtain an all-or-nothing solution if we scale the capacities C in advance as 
summarized in the following corollary. 

Corollary 2. Assume min e c e / f3 > max^ A e un. Run the GVOP algorithm with 
scaled capacities C / (3. The solution Y is an all-or-nothing 0(f3) -competitive so- 
lution. 



4.3 A Reduction of Requests with Durations 

We now add durations to each request. This means each request rj is characterized, 
in addition, by a duration interval Tj = [t*p , t^], where rj arrives in time t*p and 

ends in time Requests appear with increasing arrival times, i.e., < 
For example, the capacity constraints in virtual circuits now require that, in each 
time unit, the bandwidth reserved along each edge e is at most c e . The benefit 
obtained by serving request rj is bj ■ \Tj\, where \Tj \ = t^p — t^p . We now present 
a reduction to the general framework. 

Let r(j, t) denote a 0-1 square diagonal matrix of dimensionality J2i<j 
The diagonal entry corresponding to fa g equals one if and only if request is 
active in time t, i.e., r(j,t)ug\ug\ = 1 iff t G T; L . The capacity constraints are 
now formulated by 

Vt : Aj ■ r(j, t) ■ Yj < C. 

Since r(j, t) is a diagonal 0-1 matrix, it follows that each entry in A(j, t) = 
Aj ■ r(j, t) is either zero or equals the corresponding entry in Aj . Thus, the as- 
sumption that maxj^ A e ^^ < c e and that minjj A e ^ ^ E {0} U [1, oo) still hold. 
This implies that durations of requests simply increase the number of capacity con- 
straints; instead of Aj - Yj < C, we have a set of ./V constraints for every time unit. 
Let T max denote maxj Tj. Let Aj denote the N ■ (t^ + T max ) x J2i<j matrix 
obtained by "concatenating" A(j, 1), . . . , A{j, tf Y ), A(j, tf } + T max ). The 

new capacity constraint is simply Aj Yj < C. 

Fortunately, this unbounded increase in the number of capacity constraints has 
limited implications. All we need is a bound on the "weight" of each column of Aj. 
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Consider a column (i, £) of Aj. The entries of this column are zeros in A(j, t') for 
t' £Ti. It follows that the weight of column (i, I) in Aj equals |T$ | times the weight 

of column (i,£) in A(i, if ). This implies that the competitive ratio increases to 
(2, j3') -competitiveness, where /3' = log 2 (l+3-T max -(maxj^w(j,£))-(maxj bj)). 

Theorem 3. The GVOP algorithm, when applied to the reduction of online packing 
with durations, is a (2, (5')-competitive online algorithm. 

Remark 2. Theorem \3\ can be extended to competitiveness in time windows 
This means that we can extend the competitiveness with respect to time intervals 
[0, t] to any time window [tx, t%]. 

Remark 3. The reduction of requests with durations to the online packing frame- 
work also allows requests with split intervals (i.e., a union of intervals). The dura- 
tion of a request with a split interval is the sum of the lengths of the intervals in the 
split interval. 

Remark 4. In the application of circuit switching, when requests have durations, it 
is reasonable to charge the request "per hit" . This means that bj/ (dj ■ \TA) should 
be within the range of prices charged per bit. In fact, the framework allows for 
varying bit costs as a function of the time (e.g., bandwidth is more expense during 
peak hours). See also for a discussion of benefit scenarios. 

4.4 Approximate Oracles 

The GVOP algorithm relies on a VNet embedding "oracle" which computes 
resource-efficient realizations of the VNets. In general, the virtual network em- 
bedding problem is computationally hard, and thus Step[T]could be NP-hard (e.g., 
a min-cost Steiner tree). Such a solution is useless in practice and hence, we ex- 
tend our framework to allow for approximation algorithms yielding efficient, ap- 
proximate embeddings. Interestingly, we can show that suboptimal embeddings do 
not yield a large increase of the competitive ratio as long as the suboptimality is 
bounded. 

Concretely, consider a p-approximation ratio of the embedding oracle, i.e., 
l(j,£) < P ■ min{7(j',£) : f j;i G Aj}. The GVOP algorithm with a p- 
approximate oracle requires two modifications: (i) Change the condition in Step [2] 
to j (J, I) < bj ■ p. (ii) Change Step|2b]to Zj bj ■ p — t) / p. 

The following theorem summarizes the effect of a p-approximate oracle on the 
competitiveness of the GVOP algorithm. 

Theorem 4. Let (3 p = p ■ log 2 (l + 3 • p ■ (maxj^ w(j, £)) • ( maxj bj)). Under 
the same assumptions of Theorem\2\ the GVOP algorithm is a (2, f3 ' ^-competitive 
online all-or-nothing packing algorithm if the oracle is p-approximate. 
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5 Application to VNet Service Models 



In this section we show how the framework for online packing can be applied to 
online VNet embeddings. The key issue that needs to be addressed is the oracles 
in Line[Uof the GVOP algorithm. 

We consider the three traffic models: customer-pipe, hose and aggregate 
ingress. We also consider three routing models: multipath, single path and tree 
routing. 

Recall that (3 in Theorem [2] is the factor by which the GVOP algorithm augments 
resources. Recall that /3' is the resource augmentation if VNet requests have dura- 
tions. The following corollary summarizes our main result as stated in Theorem [T] 
The following corollary states the values of f3 and f3' when applying Theorems |2] 
and|3]to the cases described below. 

Corollary 3. The values of f3 and f3' in Theorems\2\and\3\are (3 = 0(log(|i^| • 
(max e c e ) • (max.,- bj ) ) ) and (3' = O (log ( | T max | • | E \ ■ (max e c e ) • (max.,- bj ) ) ) for any 
sequence of VNet requests from the following types: (i) customer pipe model with 
multipath routing, ( ii) hose model with multipath routing, single path routing or 
tree routing 1 , or (Hi) aggregate ingress model with multipath routing, single path 
routing, or tree routing. 

Remark 5. Our framework can handle heterogeneous VNet requests, i.e., requests 
from any of the customer service models and routing models included in Corol- 
lary\3\ Each time a request arrives, the corresponding oracle procedure is invoked, 
without disturbing existing requests. This implies that the network resources are 
shared between requests of all types. 

5.1 Proof of Corollary |3] 

The proof deals with each traffic model separately. 

Customer Pipe Model In multipath routing, an embedding of a request is a mul- 
ticommodity flow. The flow along each edge equals the bandwidth reservation 
needed to support the request. This means that, for each request rj , the set of valid 
embeddings Aj of rj consists of all the multicommodity flows specified by the 
traffic matrix and the edge capacities. For a multicommodity flow / € Aj, the 
entry A e j equals the flow /(e). The oracle needs to compute a min-cost multi- 
commodity flow in Aj, where a cost of a unit flow along an edge e equals x e . A 

'The oracle in this case does not run in polynomial time. 
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min-cost multicommodity flow can be computed by solving a linear program or by 
a using a PTAS 1551 . 

A technical issue that needs to be addressed is that the flow along an edge may 
be positive yet smaller than one, thus violating the requirement in Theorem [2] A 
key observation is that the requirement in Theorem |2] can be relaxed to A e un € 
{0} U [jp,oo). This only affects the augmentation by a constant factor. Thus, 
one can deal with this issue by peeling off such a flow, and rerouting it along 
other paths. On the other hand, the resulting oracle in this case may violate edge 
capacities. We refer the reader to [ 19 ] where an extension of GVOP that deals with 
such an oracle is discussed in detail. 

Hose Model In multipath routing, an embedding is a reservation u of bandwidths 
to edges so that every allowed traffic can be routed as a multicommodity flow. An 
entry A e u tU \ equals the bandwidth u e reserved in e for the embedding of request 
r,j. In HH, a linear programming based poly-time algorithm is presented for a 
min-cost reservation in the hose model. As in the case of the customer pipe model, 
positive bandwidth reservations of an edge might be smaller than 1. The oracle 
in this case executes the algorithm in lfT8l to obtain an optimal reservation. This 
reservation is modified so that (1) the minimum positive bandwidth reservation is 
f2(l/m 2 ), (2) edge capacities are violated at most by a constant factor, and (3) the 
cost of the embedding is at most doubled. We refer the reader to [ 19 ] where an 
extension of GVOP that deals with such an oracle is discussed. 

An efficient approximate oracle for tree routing in the hose model is an open 
problem |29l . We elaborate on a non-polynomial oracle that focuses on the online 
aspects of the problem. A Steiner tree T is feasible if and only if the bandwidth 
u(e) reserved for each edge e equals the maximum traffic that may traverse e and 
u(e) < c e . Indeed, let A e U B e denote a partitioning of the terminals Uj induced 
by the deletion of the edge e from T. The maximum traffic along e by request rj 
equals 

min < ^2 b ou t(u), ^ b in (v) > + min I n b in (u), ^ b out {v) 

Ku£A e v£B e ) K. u£A e v£B e 

The non-polynomial oracle simply returns a min-cost feasible Steiner tree that 
spans the set of terminals Uj. If no feasible Steiner tree exists, the request is re- 
jected. A similar non-polynomial oracle exists for single path routing. 

Aggregate Ingress Model An embedding in the aggregate ingress model is also 
a reservation of bandwidths so that every allowed traffic can be routed. In the 
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multipath routing model, a min-cost embedding can be obtained by a variation 
of the algorithm presented in |[T8l combined with the a modification that avoids 
small reservations as discussed previously for the customer pipe and hose models. 
Dealing with small reservations is done similarly to the way described in the hose 
model case. 

A min-cost single path routing embedding in the aggregate ingress model is 
always a tree with uniform bandwidth reservation that equals the ingress amount. 
Thus, the routing models of single paths and trees coincide in this case. This 
implies that a min-cost tree embedding is simply a min-cost Steiner tree. The 
oracle in this case proceeds as follows. (1) Delete all edges the capacity of which 
is less than the ingress amount. If this deletion disconnects the terminals in Uj, 
then reject the request. (2) Compute a min-cost Steiner tree over the remaining 
edges (see [321 and lfT3l for the best approximation to date). 

5.2 Router Loads 

So far we have focused on the load incurred over the edges, i.e., the flow (e.g., data 
rate) along an edge is bounded by the edge capacity (e.g., available bandwidth). 
In this section we also view the nodes of the network as resources. We model the 
load incurred over the nodes by the rate of the packets that traverse a node. Thus, 
a request is characterized, in addition, by the so-called packet rate. 

In this setting, each node (router) v has a computational capacity c v that spec- 
ifies the maximum rate of packets that node v can process. The justification for 
modeling the load over a node in this way is that a router must inspect each packet. 
The capacity constraint of a node v simply states that the sum of the packet rates 
along edges incident to v must be bounded by c v . 

For simplicity, we consider the aggregate ingress model with tree routing. A re- 
quest rj has an additional parameter prj that specifies the aggregate ingress packet 
rate, i.e., prj is an upper bound on the sum of the packet rates of all ingress traffic 
for request rj . 

Applying our framework requires to add a row in A to each node (in addition to 
a row per edge). An entry A V:U equals prj if the reservation u of capacities assigns 
a positive capacity to an edge incident to v, and zero otherwise. The oracle now 
needs to compute a node- weighted Steiner tree ll27l . The approximation ratio for 
this problem is 0(log kj), where kj denotes the number of terminals in request Tj. 

The following corollary summarizes the values of p and j3 p when applying 
Theorem|4]to router loads. One can extend also Theorem|3]in a similar fashion. 

Corollary 4. In the aggregate ingress model with tree routing, p = 
0(log max.,- kj) and j3 p = 0(p ■ log(p • (\E\ • (max e c e ) + |V| • (max„Cy)) • 
(max.,- bj))). 
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6 Discussion 



This article presents a unified algorithm for online embeddings of VNets. Each 
VNet request consists of endpoints, quality-of-service constraints, and a routing 
model. The algorithm can handle VNets requests from different models (e.g., the 
customer-pipe, hose, and aggregate-ingress models), and each request may allow 
a different routing model (e.g., multipath, single-path, and tree -routing). Since the 
problem we address is a generalization of online circuit switching Q, it follows 
that the lower bounds apply to our case as well. Namely, the competitive ratio of 
any online algorithm is J7(log(n • T max )), where n denotes the number of nodes 
and T max is the maximal duration. 

In the context of hose model with tree routing, we emphasize that finding a 
feasible Steiner tree is NP-hard to approximate within a constant factor [29]. One 
can relax the feasibility requirement, and compute a Steiner tree that violates the 
capacities of G by a factor of \i > 1 while satisfying the ingress/egress demand of 
every terminal. The oracle in this case will be bi-criteria. Bi-criteria oracles can 
be incorporated in the primal-dual scheme as shown in |fl9ll . To our knowledge, 
the question whether there is a polynomial time bi-criteria algorithm for min-cost 
Steiner trees in the hose model is open. 
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